TIME DURATION OF THE PARAMETRIC X-RAY RADIATION 
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Time evolution of the parametric X-Ray radiation, produced by a relativistic charged particle 
passing through a crystal, is studied. The most attention is given to the cases when the radiation 
lasts much longer (tpxR ~ 0.1 ns) than the the time t p of the particle flight through the crystal (t p ~ 
1 ps). It is shown that such long duration of the radiation makes possible the detailed experimental 
investigation of the complicated time structure of the parametric X-ray pulses, generated by electron 
bunches, which are available with modern acceleration facilities. 
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The investigation of stop-actions in the femtosecond time scale is one of the most important part of research programs 
on modern accelerators, used for X-ray lansing (FLASH [l[, XFEL 0, LCLS 0, SPring-8 g, SwissFEL [|). Being 
able to produce femtosecond electron bunches the accelerators can provide researchers with important information 
about X-ray radiation generated by charged relativistic particles moving in crystals [6|-|8| . 

It is well known the refraction of the electromagnetic field associated with the electron passing through matter 
with a uniform velocity originates the Vavilov-Cherenkov and transition radiations [9(. The Cherenkov emission of 
photons by a charged particle occurs whenever the index of refraction n(u>) > 1 (ui is the photon frequency). For 
X-ray frequencies being higher than typical atom's ones the refractive index has a universal form [9] 



where ojl is the Langmuir frequency. As a result, n(uS) > 1 and the Vavilov-Cherenkov effect is seemed to be absent 
in the X-ray range. 

But as it was first shown in [ill EH the diffraction of photons emitted by a relativistic charged particle in crystals 
leads to an abrupt change in the refractive index of X-ray quanta and, as a result, gives rise to the the spontaneous 
and induced quasi-cherenkov radiation called the parametric X-ray radiation (PXR). It also leads to a considerable 
£S| • change in spectral-angular characteristics of transition and bremsstrahlung radiation at both small and large angles 
\ to the direction of particle motion [jj], [l3| • 

In accordance with [ill E3 the diffraction of virtual photons in a crystal can be described by a set of refractive 
t-H , indices n^u), k), some of which may appear to be greater than unity. Here k is the refraction index of a crystal for the 
X-rays with the photon wave vector k. Particularly, in the case of the two wave diffraction, the refractive indices are 
ni(u>, k) > 1 and ri2(w, k) < 1, and accordingly, the two types of waves propagate in a crystal: a fast wave (n2 < 1) 
and a slow one (nj > 1). For a slow wave the Cherenkov condition can be fulfilled, but not for a fast wave. The latter 
wave is emitted at the vacuum-crystal boundary (or due to bremsstrahlung radiation at multiple scattering). 

At present there are lots of theoretical and experimental works devoted to the PXR being first registrcd experimcn- 
taly in 1985 [ll|. Moreover the investigation of the tunable sources of radiation based on the parametric radiation 
are in progress. 

It should be noted, however, that until now, theoretical and experimental analysis of radiation produced by a 
relativistic particle passing through a crystal has focused on spectral-angular characteristics of radiation. Nevertheless, 
it was shown in [6|-|8j that because of diffraction, emitted photons have the group velocity v v gr appreciably smaller 
than the velocity v of a relativistic particle. As a result, the situation is possible in which radiation from the crystal 
still continues after the particle has passed through it p-Q- This enables studying time evolution of the process of 
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photon radiation produced during particle transmission through a natural or artificial (photonic) crystal, or during 
the particle flight along the surface of such crystals In the present paper the formulas are derived, which 

describe the time evolution of the parametric X-ray radiation produced by a relativistic particle moving in a crystal. 
It is shown that the conditions are realizable under which the PXR lasts considerably longer than the time t p of the 
particle flight through the crystal. 

Time resolution of modern X-ray detectors used in FEL technology lies in the picosecond range [HI, H3| • Conse- 
quently time dependence of the parametric radiation can be investigated experimentaly. 

The paper is organized as follows. In Sec. 1 formula for the group velocity of photons in crystals is derived in 
the two wave approximation. It is shown that the group velocity can be much smaller than the speed of light in the 
vacuum. Long time reflection of the X-ray pulses from crystals demonstrates this (Sec. 2). Equations, describing 
the time evolution of the PXR, are derived in sec. 3. The dependence of X-ray pulses on geometrical parameters are 
analysed. 



I. GROUP VELOCITY OF THE ELECTROMAGNETIC WAVES IN CRYSTALS 

Let us consider the pulse of electromagnetic radiation passing through the medium with the index of refraction 
n(cj). The group velocity of the wave packet is as follows: 
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where c is the speed of light, uj is the quantum frequency. 

2 

In the X-ray range ( ~tens of keV) the index of refraction has the universal form n(u>) = 1— ir\ , lul is the Langmuir 
frequency. Additionally, n — 1 ~ 10~ 6 <C L Substituting n(uj) into one can obtain that v gr ~ c ^1 — ^ji^j - It is 

clear that the group velocity is close to the speed of light. Therefore the time delay of the wave packet in a medium 
is much shorter than the time needed for passing the path equal to the target thickness L in the vacuum. 

AT = ~ (3) 

Vgr C C U) A C 

To consider the pulse diffraction in a crystal, one should solve Maxwell's equations that describe a pulse passing 
through a crystal. Maxwell's equations are linear, therefore it is convenient to use the Fourier transform in time and 
to rewrite these equations as functions of frequency: 



- u 2 -. _ 

-curl curl E^(r,uj) + —E^(f,uj) 



+ Xii(?,to)E i J?,u)=0, (4) 



where Xjj(r, w) is the spatially periodic tensor of susceptibility; i,j — 1,2,3 repeated indices imply summation. 

Making the Fourier transformation of these equations in coordinate variables, one can derive a set of equations 
associating the incident and diffracted waves. When two strong waves are excited under diffraction (the so-called 
two-beam diffraction case), the following set of equations for determining the wave amplitudes can be obtained: 

(S-l-Xo)^-P.X-^=0 

(5) 

(S-1-Xo)4 t -c 5 x^ = 

Here k is the wave vector of the incident wave, fop = k + f , f is the reciprocal lattice vector; xoi Xf are the Fourier 
components of the crystal susceptibility: 



X 



(r) = ^x?exp(zfr) (6) 



C s = e s e/, e s (e- s ) are the unit polarization vectors of the incident and diffracted waves, respectively. 

The solvability condition for the linear system ([5]) leads to a dispersion equation that determines the possible wave 
vectors A: in a crystal [U E3- ^ is convenient to present these wave vectors as: 
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where = 1,2; N is the unit vector of a normal to the entrance surface of the crystal, which is directed into the 
crystal, 

4 1,2) =\[0- + A)xo - Pia B ] ± \ {[(1 + A)xo - /3i«s - 2 Xo ] 2 + 4^C| X ^X-r} 1/2 , (8) 

as = (2kf + T 2 )k~ 2 is the off-Bragg parameter («g = when the Bragg condition of diffraction is exactly fulfilled), 

7o = Hy-N, n 7 = \, j8 = 71 = n~ fT -N, n 7r = , C s = 1 (ex-polarization), = cos(26>_b) (7r-polarization, 

6*b — Bragg angle). 

The general solution of equations flU, © inside a crystal is: 
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§(r) = k exp^fc^f) + e T s A TM exp^A^f) (9) 



Associating these solutions with the solutions of Maxwell's equation for the vacuum area one can find the explicit 
expression for Ei(f) throughout the space. It is possible to discriminate several types of diffraction geometries, namely, 
the Laue and the Bragg schemes, which are most well-known [l7i j . 

In the case of two-wave dynamical diffraction, the crystal can be described by two effective indices of refraction 



4 1,2) = \ {X0(1 + Pi) ~ Pl*B ± - K) + PlU B ) 2 + A^CsXrX-r} ■ (10) 

The diffraction is significant in the narrow range near the Bragg frequency, therefore xo and Xr can be considered 
as constants and the dependence on ui should be taken into account for as = k ^ = —-p-^(uj — lob), where 

k = — ; t is the reciprocal lattice vector which characterizes the set of planes where the diffraction occurs; Bragg 
frequency is determined by the condition as = 0, i. e. r 2 + 2 k b~^ = and ljb = — %=f c ^t ■ 

From ©, © one can obtain 
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(ii) 
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From ([TT|) it follows, that the complicated character of the wave field in a crystal leads to both the negative and 
the positive values of the group velocity. Moreover when 0i is negative the radicand in (|11[) approaches zero (Bragg 
reflection threshold) and v gr — > . The performed analysis allows one to conclude that the center of the X-ray pulse 
in a crystal can undergo a significant delay AT 3> — available for experimental investigation. Thus, when /3j = 10 2 , 
L = 0.3 mm and L/c ~ 1 ps, the delay time can be estimated as AT w ~ 0.1 ns. 

II. WAVE-PACKET REFLECTION 



Let us study now the time dependence of intensity of the wave-packet after passing through a crystal more datail. 
Assuming that B{ui) is the reflection or transmission amplitude coefficients of a crystal, one can obtain the following 
expression for the pulse form 
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Fig. 1. Group velocity of the slow wave for the Bragg diffraction case in LiF crystal (/3i = —45). 



E{t) = ^- J B(uj)E (uj)e- iuJt duj = J B(t - t')E (t')dt' . (12) 

where Eq(ui) is the amplitude of the electromagnetic wave incident on a crystal. 

In accordance with the general theory, for the Bragg geometry, the amplitude of the diffraction-reflected wave for 
the crystal width much greater than the absorbtion length can be written as [l7| : 

B.(u) = (13) 
{xo(l + IAD - \Pi\a B - V(Xo(l - 1/3,1) - Iftl a B f - 4 \^ \ C sX rX-r } 

In the absence of resonance scattering, the parameters xo and \±t can be considered as constants and frequency 

2 

dependence is defined by the term as = — -p— (^ -wb). 
So, B s (t) can be found from 

B ^ = -4^7 < 14 > 

x j {xo(l + - \0i\a B - V(xo(l - 1/3x1) - |/3| - 4 |AI C.XrX-r} e^dw. 



The Fourier transform of the first term results in S(t) and we can neglect it because the delay is described by the 
second term. The second term can be calculated by the methods of the theory of function of complex argument: 



B s (t) = --L !/3l | I J ^st) e _ i{UB+AuB)tm (15) 

4Xr k%LO B t 



or 



Bs (t) = Jifat) e . i{UB+AuB)te , t y (lg) 

2 a s t 



where 



2^C s XtX-tUb . Xo(l+ |ySi|)wsA;| 

fl S = i^B = ,o | 2 
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Since xo and Xt are complex, both a s and Aug have real and imaginary parts. According to (|14j) - (|15[) . in the case 
of Bragg reflection of a short pulse (the pulse frequency bandwidth ^> frequency bandwidth of the total reflection 
range) both the instantly reflected pulse and the pulse with amplitude undergoing damped beatings appear. Beatings 
period increases with grows and Xt decrease. Pulse intensity can be written as 



i s (t) ~ \B s (t)\ 2 = IJ y- 



Ji(a s t) | 



at 



-2ImAw B t 



0(t). (17) 



It is evident that the reflected pulse intensity depends on the orientation of photon polarization vector e s and 
undergoes the damping oscillations on time. 

Let us evaluate the effect. Characteristic values are ImAuiB ~ Imxo^B and Ima ~ . For 10 keV for the crystal 

of Si Imxo = 1,6- 1(T 7 , for LiH Imxo = 7, 6 • 10" 11 , Imxr = 7 • 10" 11 , for LiF Imxo ~ 10~ 8 . Consequently, the 
characteristic time r for the exponent decay in (fTT5)) can be estimated as follows (wg = 10 19 ): 

for Si the characteristic time r ~ 10~ 12 sec, for LiF the characteristic time r ~ 10~ 10 sec, for LiH the characteristic 
time r ~ 10 -9 sec!! 

The reflected pulse also undergoes oscillations, the period of which increases with growing f3i and decreasing Rexr- 
This period can be estimated for (3 = 10 2 and Rex T ~ 10~ 6 asT~ 10 -12 sec (for Si, LiH, LiF). 

When the resolving time of the detecting equipment is greater than the oscillation period, the expression (| 16[) should 
be averaged over the period of oscillations and the delay law (JTHJ) has the power function form: 

i s (t) - rV^'W. 

In the case of multi-wave diffraction, the time delay for the photon exit from the crystal will be even more. 



III. PARAMETRIC X-RAY RADIATION 



So the time t p u — that photons spend in a crystal can be longer than the flight time t p = ^ of a relativistic 
particle in a crystal. Hence, the emission of diffraction-related PXR (diffracted radiation of an oscillator, surface 
parametric radiation and others) produced by a relativistic particle will continue after the particle has left the crystal. 
Under diffraction conditions, the crystal acts as a high-quality resonator. 
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Fig. 3. Schemes of diffraction in the Bragg geometeries: diffracted maximum (left), forward maximum (right). 



The intensity I(t) of radiation produced by a particle which has passed through a crystal can be found with known 
intensity of the electric field E(f,t)) (magnetic field H(f,t)) of the electromagnetic wave, which is produced by this 
particle @, 



I(t) = -\E(r,t)\ 2 r 2 , 

where r is the distance from the crystal, which is assumed to be larger than the crystal size. 
The field E(f, t) can be presented as an expansion in a Fourier series 



(18) 



(19) 



At a long distance from the crystal, the Fourier component E(f, ui) can be written as follows 

Ei(f,u) = — %Y, e ° 1 4 _) *V")M")dV. (20) 



where i = 1,2,3 (and correspond(s) to the coordinate axes x, y, z), e| is the z-component of the wave polarization 



~ > (—)s 

vector e 8 : s = 1, 2: El is the solution of Maxwell's equations, 

k 



j(f,cj)= / j(r,t)e iut dt 



(21) 



j(f,u>) = Qv(t)S(f— r(t)) is the is the current density of the particle with charge Q, r(t) is the particle coordinate at 
time t. 

The explicit form of the expressions E\ describing diffraction of the electromagnetic wave in a crystal in the 
Laue and Bragg cases is given in It should be noted there is the close connection between Si ' s and E^ s , 

describing scattering of a plane wave with a wave vector k = on a target and the asymptotic of a diverging spherical 
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wave: E<f )s = E ( ^ s * p^4l6| . This circumstance permits us to use the results obtained in diffraction theory to find 
explicit expressions for e\ 

k 

Now let us look at the expression for the amplitude A(ui) of the emitted wave more attentively: 

A' n (u,) = % I Ef )s (r'.^)/(r'..') ( /-V'. 122) 



Using (|2ip . (1^1) can be recast as follows 



AUu) = ^ I E ( f )s ' (7,u)Qv(t)6(i* -r(t))e iut dtd 3 r' 



^ f Et )s '{r{t),u)v{t)e wt dt (23) 

According to f|23[) . the radiation amplitude is determined by the field E~ ^ taken at point r(t) of particle location 
at time t and integrated over the time of particle motion. 

From (fH?|) , , and (1^1) it follows that the expression for the electromagnetic wave emitted by the particle passing 
through the crystal (natural or photonic) can be presented in a form: 

E^t) = ^E e "/ At^e-^-Vdu, (24) 

i.e., = i£eMc(*-§)- 

S 

Consequently the expression for the X-ray intensity can be written as 

mt) = £\^At(t-l)f. (25) 

S 

If we are interested in the intensity of emitted photons with definite polarization then 

W,t) = ^\<?Al(t- r ~)\ 2 . (26) 

From ( (|24)l ) follows that the time dependence of the form of the pulse I(f,t)(E(f,t)) of radiation generated by a 
particle passing through the crystal is determined by the dependence of the radiation amplitude Ai(u) on frequency. 
According to (|23|) to find explicite expressions for the radiation amplitudes Ai(ui) it is neccesary to know wave 

functions E~ ^ s , which have been obtained in 



A. Laue case 



Let us consider Laue case (fig. 2). In this case electromagnetic waves emitted by a particle in the forward and 
diffracted directions leave the crystal through the same surface. By matching the solutions of Maxwell's equations on 
the crystal surfaces one can obtain the following expression for the amplitude in Laue case for forward directions [l5[ 



= E & 



i — e us L 
e M 



li=l,2 



uj — kv uj — (k + k^ s N)v 



(27) 



and for diffracted directions 



E e. 

1 



u> — k T v uj — (k T + K^gN)! 



(28) 
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Fig. 4. The PXR in diffracted direction in the Bragg geometry in LiF crystal for electron energies 10 GeV (right), 1 GeV 
(left): Pi — —45 (solid line), fii = —30 (dashed line) (diffraction plate (5,5,5)). 



where 



k T = k + 

aj __ 2e 2(i) 
sn 

Cl(2)s 



Xo 



^- + 2( £2s - £ls ) 
_ C s x-f 
2(e 2s - eis)' 



(29) 



B. Bragg case 

Now let us consider the PXR in Bragg case (fig. 3). In this case, side by side with electro-magnetic wave emitted in 
the forward direction, an electromagnetic wave emitted by a charged particle in the diffracted direction and leaving 
the crystal through the surface of particle entrance can be observed. By matching the solutions of Maxwell's equations 
on the crystal surfaces we can get the formulas for the amplitude in Bragg diffraction case for forward directions 



U! — kv LU — (k + K^sN)!! 



J(u}—(k+K, fis N)u)- 



and diffracted directions 



H=l,2 



i(ui— (k T + n lis N)u)- 



whcrc 



7l(2)a 



7l(2)s 



iq — k T v uj — (k T + k^ s N)v 

2e 2 (i) s - Xo 



2e 2 (i) s - Xo - (2ci( 2 ) s - Xo) exp (^(e 1(2)s - e 2 (i )s )£) ' 

-PiC s Xf 

2e 2 (i) s - Xo - (2e 1(2 ) s - Xo) exp (^-(ei( 2 ) s - ^2{i)s) L ) 



(30) 



(31) 



(32) 



With the help of expressions (|2"T|) — (|3"2")l time properties of the PXR can be investigated. 

In the subsequent analysis it should be taken into account the parameter |ob| ~ — ^>b\ rising the PXR intensity 
decreases rapidly, i. e. the radiation is concentrated near Bragg frequency u>b- This gives us an opportunity to 
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represent the rate of photon emission in the following form 

d 2 N _ I(f,t) 
dtdfl Itwb 



(33) 



The time evolution of the PXR in diffracted direction is introduced on fig. 4. Due to the small values of the group 
velocity in the antisymmetric diffraction scheme the time of emission lies in the nanosecond range. Moreover the 
radiation pulse has complicated structure: it undergoes damping oscillations. The period and the height of radiation 
peaks increase with f3\. 

It should be noted according to time duration of the quasi cherenkov radiation produced by relativistic particles 
in optical and microwave ranges in the case XOtXt 1 describe by the formulas analogous to the above. 

It should be mentioned that the above discussed phenomena can not be described in the framework of the PXR 
theory using the first Born approximation (though there were a lot of attempts to apply it for consideration of the 
PXR spectral-angular properties). 

For observation of oscillations, one should either register the moment of particle entrance into the crystal or use a 
short bunch of particles with duration much shorter than the oscillation period. In the X-rayrange, such situation 
is typical of electron bunches, which are applied for creating X-ray FELs (FLASH XFEL LCLS Q, SPring-8 
[H, SwissFEL Q). (The bunch duration in such FELs is tens-hundreds of femtoseconds). If the bunch duration is 
large in comparison with the duration of the radiation pulse or the time of the electron entrance into the crystal is 
not registered, which occurs in a conventional experimental arrangement, then the intensity I(t) should be integrated 
over longer observation time intervals. As a result, we, in fact, obtain the expression integrated over all frequencies, 
i.e., an ordinary stationary angular distribution of radiation. If the response time of the devices detecting tjj (or the 
flight time of the particle in a crystal t p , or the bunch duration) is comparable with the oscillation period, then I{t) 
should be integrated over the interval to- In this case oscillations will disappear, but we will observe the decrease in 
the intensity of radiation from the crystal. 



Conclusions 



The formulas which describe the time evolution of radiation produced by a relativistic particle moving in a crystal 
are derived. It is shown that the conditions are realizable under which the parametric X-ray radiation lasts much 
longer than the time t p of the particle flight through the crystal. It is shown that such long duration of the radiation 
makes possible the detailed experimental investigation of the complicated time structure of the parametric X-ray 
pulses, generated by femtosecond electron bunches, which are available with modern acceleration facilities. 



[1] flash.desy.de/ 

[2] xfel.desy.de/ 

[3] lcls.slac.stanford.edu/ 

[4] www-xfel.spring8.or.jp/ 

[5] fel.web.psi.ch/ 

[6] V. G. Baryshevsky // Izvestia AN BSSR. Ser. Phiz.-mat. 1989. No. 5. P. 109. 

[7] V. G. Baryshevsky // arXiv:9906022vl ( 1999). 

[8] V. G. Baryshevsky, A. A. Gurinovich // |arXiv:1101.4162H (2011). 

[9] L. D. Landau, E. M. Lifshitz, The Electrodynemics of Continues Media in: L.D. Landau, E.M. Lifshitz Course of Theoretical 
Physics Vol. 7 (Pergamon Press, 4ed., 1975). 
[10] L. D. Landau, E. M. Lifshitz, The Classical Theory of Fields in: L.D. Landau, E.M. Lifshitz Course of Theoretical Physics 

Vol. 2 (Pergamon Press, 4ed., 1975). 
[11] V. G. Baryshevsky // DAN BSSR. 1971. Vol. 15. P. 306. 

[12] V. G. Baryshevsky, I. D. Feranchuk // Zh. Exp. Teor. Fiz. 1971. Vol. 63. P. 944. 

[13] V. G. Baryshevsky , O. M. Lugovskaya // Physics of Atomic Nuclei. 2003. Vol. 66. P. 409. 

[14] V.G. Baryshevsky Channeling, Radiation and Reactions in Crystals at High Energy (Bel. State Univers., Minsk, 1982) 
[15] V. G. Baryshevsky // NIM B. 1997. Vol. 122. P. 13-18. 

[16] V. G. Baryshevsky, I. D. Feranchuk, A. P. Ulyanenkov, Parametric X-Ray Radiation in Crystals: Theory, Experiment and 

Applications (Series: Springer Tracts in Modern Physics, Vol. 213 2005). 
[17] Z. G. Pinsker Dynamical scattering of X-rays in crystals (Springer, Berlin, 1988) 

[18] Chang Shih-Lin, Multiple Diffraction of X-Rays in Crystals ( Springer- Verlag Berlin Heidelberg New- York Tokyo, 1984) 

[19] Z. Chang et al. (1996), App. Phys. Lett. 69, 133-135 

[20] P. Gallant et al. (2000), Review of Scientific Instruments 71, 3627-3633. 



